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We analyze a microscopic mechanism behind coexistence of a heavy Fermi liquid and geometric 
frustration in Kondo lattices. We consider a geometrically frustrated periodic Anderson model and 
demonstrate how orbital fluctuations lead to a Kondo-screened phase in the limit of extreme strong 
frustration when only local singlet states participate in the low-energy physics. We also propose a 
setup to realize and study this exotic state with SU (3)-symmetric alkaline-earth cold atoms. 
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Introduction. Geometric lattice frustration plays a 
crucial role in Mott insulators [1] where it usually sup¬ 
presses long-range magnetism by enhancing the number 
of competing magnetic ground states. At zero tempera¬ 
ture, this degeneracy may be relieved in favor of a quan¬ 
tum non-magnetic phase such as a spin liquid or valence 
bond ordering [2]. On the contrary, lattice topology in 
most metals is less important due to long-range magnetic 
interactions mediated by the itinerant electrons and small 
static magnetic moments. 

The situation is different in cases when magnetic and 
itinerant behaviors originate from physically distinct de¬ 
grees of freedom [3]. For example, in heavy-fermion 
(HF) metals [4, 5] magnetic moments arise from localized 
4/ or 5/-electrons, while conduction electrons typically 
reside in extended atomic s-orbitals. Low-temperature 
properties of such systems are driven by several oppos¬ 
ing quantum many-body effects: (i) Kondo screening, 
i.e. formation of singlets between local moments and 
itinerant electrons that gives rise to “heavy” quasiparti¬ 
cle states with delocalized /-electrons; (ii) local-moment 
long-range magnetism; and (iii) non-magnetic states due 
to lattice frustration that involve singlets only among 
local spins. Geometrically frustrated /-electron com¬ 
pounds [or Kondo lattices (KLs)] such as Yb 2 Pt 2 Pb [6] 
received much attention in the recent years [7-10]. 

The magnetism, Kondo effect, and geometric frustra¬ 
tion compete because they involve same local electrons 
which can not simultaneously form singlets with each 
other and the conduction band. This observation is at 
the heart of the recently proposed generic phase diagram 
of HF materials [7] that allows their classihcation ac¬ 
cording to the amount of quantum fluctuations of local 
magnetism [11]. Naturally, this phase diagram precludes 
Kondo screening in strongly-frustrated lattices. The an¬ 
tagonism between Kondo effect and lattice frustration 
only occurs in cases that involve pure spin degrees of free¬ 
dom. In contrast, in systems with multiple local orbitals, 
orbital fluctuations allow local spin singlets to participate 
in the Kondo screening [12, 13] together with the usual 
“spinful” states. If the singlets were due to frustration, 
the local orbital fluctuations might provide a pathway 


towards a strongly-frustrated Kondo-screened state. 

In the present Letter we argue that such phase with 
coexisting Kondo and frustration-driven local-spin sin¬ 
glets can indeed be realized. To demonstrate this, we 
consider a toy system - a periodic Anderson model on a 
triangular tube lattice (TTL) of Fig. 1(a) with frustrated 
triangular plaquettes (due to large antiferromagnetic ex¬ 
change interaction between localized electrons) in the 
Kondo regime when valence fluctuations are suppressed 
and each plaquette has a spin-singlet ground state (GS) 
with exactly two fermions. Because of different possible 
arrangements of local valence bond (VB) singlets [14], 
this GS is triply degenerate. Although local spins are 
quenched in the singlet states, orbital fluctuations [Fig. 
1(b)] allow mixing of the VB configurations by the An¬ 
derson hybridization with the conduction band, and give 
rise to a robust Kondo-screened GS with heavy quasipar¬ 
ticles and delocalized VB singlets [Fig. 1(c)]. 

This KL can be implemented using fermionic alkaline- 
earth atoms (AEAs), i.e. atoms with two outer electrons, 
in an optical lattice [see Fig. 1(d)]. AEAs prepared in the 
two lowest clock states (^S'o and ^Pq) with total angular 
momentum J = 0 show a strict decoupling of electronic 
orbital and nuclear-spin degrees of freedom, and obey an 
accurate SU{N ^21+1) (/ is the nuclear spin) sym¬ 
metry in the two-body collisions [15] which has been re¬ 
cently verified with ®^Sr [16] and ^^^Yb [17, 18]. Our key 
observation is that the local VB singlets can be encoded 
with entangled states of two AEAs [Fig. 1(e)] prepared in 
different clock configurations and three nuclear spin lev¬ 
els. The degeneracy of these states is guaranteed by the 
SU{N = 3) symmetry. The entangled atomic pairs are 
loaded in the lowest, strongly localized, band of a magic 
optical lattice whose trapping potential does not affect 
clock transitions [19], and implement the locally frustrated 
plaquettes (the optical lattice itself does not need to be 
geometrically frustrated). The conduction electrons are 
simulated by placing AEAs in a higher, itinerant band. 

At low energies, both of the above systems are de¬ 
scribed by a KL model with a peculiar SU (3) structure. 
In the metallic regime, its GS is a Fermi [in one dimen¬ 
sion (ID), Luttinger] liquid consisting of delocalized VB 
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FIG. 1. Panel (a) The triangular tube lattice. Black and blue 
circles denote local (dj^^) and itinerant electrons. Grey 

ellipses are VB singlets (2) (empty circles indicate holes). Itin¬ 
erant fermions propagate (by hopping between triangles with 
an amplitude t) in the leads with momentum k. (b) Valence 
fluctuations away from the two-electron singlet GS of a trian¬ 
gle, leading to VB flips, (c) Schematic plot of delocalization 
of the VB singlets and heavy-quasiparticle formation due to 
the Kondo screening. Shaded regions are Fermi surfaces in 
the Brillouin zone, (d) Magic optical lattice that implements 
the frustrated KL model (3) [band 1 (2) is localized (itin¬ 
erant)]. Red and blue circles are AEAs in (e) and 
(g) clock states. Grey ellipses show local e-g entangled states 
(with energy — Veg). Peg and Vgg are the inter-band exchange 
interactions (Kg S> Kg, IKsD- (®) Mapping from VB sin¬ 
glets on a triangle to lowest-energy e-g pair states for SU (3) 
AEAs. (f) Two different scattering lengths for spin symmetric 
\t} and antisymmetric js) channels. The three levels indicate 
nuclear spin states for each atom (black circle marks a popu¬ 
lated state). 


singlets (AEA pairs) screened by itinerant fermions, that 
can be viewed as a short-range resonant VB spin liquid 
[20] stabilized by the Kondo effect. 

Toy model: SU{3) KL on a TTL. Let us consider a 
periodic Anderson model on the lattice of Fig. 1(a): 


HTTL = Ekacl^^Ckac^ + V + h.c)-b (1) 


ka<7 


CdNfifci) -\- 

i ^ a^b 


H„ 


which describes a system of conduction electrons Ckaa 
with momentum k in the ath lead (a = 1...3), spin 
= {t, i}, hybridized (via an amplitude v) with lo¬ 
cal electrons diaa at each vertex a of a triangle at po¬ 
sition X^ = i. Nd = J2anfa = J2aa4aadiaa and 


Sd = I J^aa/s 4aa'^«/3d*a/3 (cT are Pauli matrices) define 
electron number and total spin of a triangle. The dis¬ 
persion Ska- = Ck — her includes a small (compared to 
other magnetic interactions) Zeeman splitting h whose 
role we explain later. The term iLmix describes mixing of 
fermions in different leads a and for now will be ignored. 

There are several energy scales associated with each 
triangle: local binding energy Cd > 0, the nearest- 
neighbor Coulomb repulsion U, “Hund” energy Jh ^ 0 
that forces the lowest total spin Sd, and an infinitely large 
on-site Coulomb repulsion preventing double occupancy 
of any vertex a. We focus on a two-electron Sd = 0 sub¬ 
space which contains a three-fold degenerate GS when 
U — \ Jh < ed < 217 -b \ Jh'- 

I®)* = 72 (2) 

where s^,f, = s^/,, = 1 when a, b and b' are different, and 
0 otherwise; jvac) is the vacuum (Nd = 0) state. These 
states are labeled by the number of an unoccupied vertex. 

We will fix = fC' and consider the strong-coupling 
regime u <C U, Jh when A^-fluctuations on each tri¬ 
angle are virtual and can be taken into account via a 
generalized Schrieffer-Wolff transformation S [21] that 
includes processes shown in Fig. 1(b). A straightforward 
calculation yields the second-order KL Hamiltonian 

Hgf ='Y^^kij4aaCkaa — 'Y^^<rbfiafib4aa^ibiy (3) 

ka<j iaab 

that describes scattering of conduction electrons by the 
local VB singlets and is defined on a non-frustrated lattice 
whose sites correspond to triangles in Fig. 1(a). The 
coupling constants are Vab = Kl(1 — Sab) + Vi\Sab with 
y± = y\\ = Oab - the Kronecker delta, and the 

valence fluctuation gap A = -\- ^Jh- The states (2) 

are described with a pseudo-fermion representation [5]: 

l«)* ^ /iaivac) (4) 

with a Hilbert space constraint J^afLfia = 1- Because 
only Sd = 0 triangle states are involved in the low-energy 
physics, interactions in preserve electron spin a 

and only change the orbital (lead) degree of freedom a. 

As a result, Eq. (3) describes a two-channel KL model 
(spin is the channel index) [22]. It is known that the 
two-channel fixed point is usually unstable w.r.t. chan¬ 
nel asymmetry [22] controlled by the Zeeman splitting h. 
Since even for small h Jh the leads may be consid¬ 
ered spin-polarized, below we omit the spin index a and 

replace Ciaa ^ = Cia and Ska — ~ ^k- 

The Hamiltonian (3) contains matrix elements con¬ 
necting all three possible local VB states and conduc¬ 
tion electron “flavors” a, and is an anisotropic (XXZ- 
like) SU(3) KL model written in terms of generators 
Ta(xi) = flafib and T^(xi) = cJjCia for local and itin¬ 
erant degrees of freedom [23]. The local SU(3) “spin” 
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operators Tj^ixi) describe orbital fluctuations in Eq. (1) 
that flip the VB singlets (2). in Eq. (3) is in¬ 
variant under C/(l) transformations fia —>■ fia and 

Cia —t that preserve the V± term. There is also 

a discrete lattice symmetry Csy = {Ca, cr«} [24] that con¬ 
tains ^ (C 3 ) rotations around the TTL axis and three 
symmetry planes dy of the triangles. 

Kondo effect-assisted VB phases. To demonstrate 
that the model (3) has a Kondo screened GS, we use 
a generalized hybridization mean-field (HME) approach 
[25] that treats the /-fermion Hilbert space constraint on 
average, jfJfiaifLfia) = 1 is the system size), and 
self-consistently compute the hybridization and SU(fi) 
“magnetization” order parameters (OPs) [26]. We as¬ 
sume that all OPs are site-independent. There are 
three hybridization amplitudes: xo = ^Y^aUiaCia), 

Xi ,2 = fi2Ci2+uj^^ faCa) with w = 

and eight magnetizations [ml] for c- [/-] fermions 
defined via (c^Cib) = + ^Sab [{fLU) = 

Yi Yb'nT'l + |<^ab] where A* are the Gell-Mann matrices, 
I = 1... 8 and is the conduction band filling. Unlike 
the real SU{2) magnetization, do not break time- 
reversal invariance but rather the above 17(1) and Csy 
symmetries. The OPs xi and X 2 are connected (up to a 
phase) by the planes (j„ from C^y. Finite mfl completely 
break C^y leading to nematic states; with 1 7 ^ 3, 8 



FIG. 2. T — 0 phase diagram of the KL model (3) with 
N = 4900 sites and electron density n‘^ = 0.8. mtot = -\- 

mlff' -\- (mg -l-rrig )^]^1^ plays a role of the total magnetization. 
The black circle at V\\/t ~ 0.76 marks the first order transition 
between Kondo-screened and normal phases. Inset: The HF 
band-structure. The band splitting at the Fermi level is ~ 
UXi, 2 - The table shows symmetries broken by different OPs 
and corresponding VB liquid phases. 


also break the above 17(1) symmetry. 

Kondo-screened states correspond to nonzero values of 
either hybridization OP Xo,i, 2 - In analogy to spin sys¬ 
tems [27, 28], we call phases with xi 7 ^ X 2 chiral [29]. 
This discussion is summarized in the table in Fig. 2. 

The phase diagram of the Hamiltonian (3) is shown 
in Fig. 2 for = —21 cos A: (1 is the nearest-neighbor 
hopping). There is a first order transition between a 
normal state with xo,i ,2 = 0, and a Kondo screened phase 
with Xi 7^ X2 7^ 0 (but Xo = 0) and non-zero TO3 §. This 
chiral nematic phase has delocalized VB singlets. The 
OPs survive only at low temperature T ^ ~ 
5 X 10“^1; for T > Ty the only finite OP is xi and the 
GS realizes a chiral metallic VB spin liquid. 

This phase is quite different from the Kondo-stabilized 
spin liquid of Ref. 30 where the resonating VBs of local 
spins are formed due to their coupling to conduction band 
and are unstable away from the Kondo regime. In our 
case the VB singlets are due to geometric frustration, and 
the Kondo screening only injects them into the Fermi sea. 

Stability of the Kondo-assisted VB liquid. The Kondo 
phase in Fig. 2 is quite robust against changes in 
the noninteracting itinerant density of states (DOS). 
To show this, we consider a model DOS that corre¬ 
sponds to a square lattice with nearest-neighbor hopping 
Ck = —2t{coskx -b cosky) (as opposed to the ID tight- 
binding dispersion used before), see inset in Fig. 3. The 
phase diagram obtained by applying the HMF approach 
to the KL (3) is presented in Fig. 3. Unlike the ID case 
in Fig. 2, the chiral VB liquid with xi 7 ^ 0, xo ,2 = 0 
and = m-I = 0 exists even at T = 0 for VI ^ 0 
and large V\\. Only mirror symmetry from C^y is broken 
by this state. With decreasing |Ul| and V|| the system 
undergoes a transition to a nematic metallic state with 
jT^c,/ ^ Q completely broken C^y. The situation is 
different for Vj_ 0. Here the only non-zero OP is xo 
and the VB liquid GS does not break any discrete sym¬ 
metry. All these Kondo-screened states become unstable 
at small |V1| and V||. 

The phase transitions in Fig. 2 and 3 are first order 
which may be an artifact of the HMF approximation. In 
general at T = 0 the emergence of nonzero OPs Xo,i ,2 
is associated with a phase transition (as opposed to a 
crossover) when fluctuations beyond HMF are taken into 
account [31]. Therefore salient features of our phase dia¬ 
grams should remain unchanged. 

Finally, we mention effects of a finite lead-mixing iJmix 
in Eq. (1). Its simplest form (compatible with Csy sym¬ 
metry of the TTL) corresponds to hopping of itinerant 
and local fermions around the triangle. This correction 
results [32] in a Zeeman-like term, proportional to the 
intra-triangle hopping, which lifts degeneracy of the lo¬ 
cal VB states (2) and can suppress the Kondo phase in 
Fig. 2 if this splitting is sufficiently large [4]. 

Implementation with ultracold AEAs. We propose 
an experimentally accessible implementation of the KL 
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FIG. 3. Generic T = 0 phase diagram of Eq. (3) [or 
(7) with Ugg = 0] with N = 3600 sites and rf = 0.8. 
All phase transitions are first order. Inside the nematic 
phase there is a “metamagnetic” transition between states 
with 7 ^ 0, JTig’g = 0 (at smaller I/ji) and rn{’^ = 0, 

mg’g 7 ^ 0 (for larger V||) which are separated by a continua¬ 
tion of the dashed line. Inset: Noninteracting itinerant DOS 
g{t) = — (e/dt)^) [K{x) is elliptic integral of the 

1st kind] used to compute the phase diagram. 


model (3) with AEAs in an optical lattice, that is free 
of the mixing described by The key idea of our 

approach is to use nuclear spin states of the atoms as 
“synthetic” frustrated plaquettes [corresponding to tri¬ 
angles in Fig. 1(a)] and construct an appropriate low- 
energy model that takes into account these local states 
as well as the itinerant degrees of freedom, and is unitar- 
ily related to the KL model (3). The GS degeneracy of 
a synthetic plaquette is guaranteed by the SU{N = 3) 
symmetry of the AEAs. 

Consider a two-band optical lattice schematically 
shown in Fig. 1(d). The lowest-energy band is local¬ 
ized and contains two AEAs per site in different clock 
states: one (GS, g) and one ^Pq (excited state, e). 
To minimize lossy e-e collisions, the higher-energy itiner¬ 
ant band is populated only by g atoms. The Hamiltonian 
of the system is [15]: 


= -tJ2icl„Cjn + h.c.) -f 
(ij) » 


^<«-l)+ (5) 


T (.^gg9in9im + ^eg cl^Cin + Veg ^\n^im9im9v 


where gin (cin) denote g (e) fermions in the localized 
band at site i and nuclear spin state n = 1, 0, 1 [h = —n, 
i.e. T = —1, 0 = 0], and create itinerant g atoms. 
There is an implicit summation over nuclear spin indices. 
The first term describes nearest-neighbor hopping with 
an amplitude t. The second sum corresponds to e-g {Veg 
and Vfg) and g-g iVgg) exchange couplings, as well as 


direct g-g interaction Ugg ^ 0 [see Fig. l(d)-(f)]. Veg 
and Vfg have the same sign, and Vgg is negative [33]. 

States of a localized e-g pair are described by the term 
H\oc{xi) = Vegglngimeln^Cin whose spectrum consists of 
a triply-degenerate GS subspace with energy —Veg'- 

\^)i ~ ~^^lnmeln9imW^^) (®) 

{sinm is the antisymmetric Levi-Givita tensor, = 1), 
and a sextet \l){ = {1/V%sin^e\^gljvac) [s[]^ was de¬ 
fined in (2)] and f)” = e\ig\f\Ya.c) with energy +14g- We 
assume that Veg is large, Veg | 14 sl> Kg’ ^ P4], ne¬ 
glect mixing of the above sectors, and project the Hamil¬ 
tonian (5) onto the subspace (6). Using the relations 

~ '^^^\9\n9'iTa\p') i — ^ral^np) ^ and 

the pseudo-fermions (4), we obtain an effective model 


^ ^kc\iCkl ^ 


kl 


vfy^p4u^p -%<«-!) 


( 7 ) 

with V = . If the states (6) are identified with VB 

singlets (2) on a triangle [Fig. 1(e)] by assigning a nuclear 
spin flavor m to each vertex, in ID is equivalent to 
(spin-polarized) in Eq. (3) with Ul = V|| = U [35] 

plus a Hubbard term, whose role as well as possible ways 
to introduce anisotropic couplings in Eq. (7) we discuss 
below. To reach a Kondo screened GS one must have 
U > 0, i.e. Vfg > -Vgg > 0 [36]. 

Discussion. Our theory highlights the fundamental 
role played by the orbital degrees of freedom in stabilizing 
a Kondo-screened phase in the presence of extreme strong 
frustration when only singlet local states participate in 
the low-energy physics, by allowing the conduction elec¬ 
trons to dynamically flip the VB singlets [see Fig. 1(b)]. 
These microscopic processes lead to delocalization of the 
local VBs and drive the formation of the VB spin liquid 
with HE quasiparticles. We illustrated this mechanism 
by studying a periodic Anderson model on a frustrated 
triangular tube, and proposed a optical lattice setup to 
realize this toy model with S'[/(3)-symmetric AEAs that 
employs their nuclear-spin degrees of freedom to imple¬ 
ment geometrically frustrated plaquettes (e.g. triangles). 

Compared to the electronic KL (3), the low-energy 
model for AEAs (7) has several peculiarities. First, there 
is the Hubbard term Ugg which below half-filling en¬ 
hances phases with non-zero SU{5) magnetization in Fig. 
3(b). However, its magnitude is effectively damped by 
the density prefactor ~ (n°)^. We checked that even 
when n° = 0.8, one needs Ugg > U to suppress the 
Kondo-screened state. Hence this term is unimportant 
for the Kondo physics. Second, the Hamiltonian 
has full SU{3) symmetry (i.e. Vj_ = V||) that originates 
from the symmetry of Eq. (5) and prohibits experimen¬ 
tal exploration of the phase diagram in Fig. 3. This 
symmetry can be broken by a weak external magnetic 
held B which to the lowest order amounts to replacing 
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V± —>■ V± + {fie — ^^g)B {fie,g are magnetic mo¬ 

ments for e and g atoms). Also, one might use alternative 
implementations of the SU{3) Kondo effect, e.g. using 
orbital degrees of freedom [37], instead of the AEAs setup 
discussed here. 

The HF phase in Figs. 2 and 3 can be detected in cold- 
atom experiments using slow quantum dynamics or time- 
of-flight measurements [38-40]. The KL model in Fq. (7) 
can be implemented beyond ID, which enables us to use 
AFAs as controlled [because of the SU{N) symmetry] 
quantum simulators for more complex frustrated Kondo 
lattices. Although the currently available isotopes ®^Sr 
and i^^Yb are believed to have negative exchange cou¬ 
plings V [16-18], we expect that our results summarized 
in Figs. 2 and 3 can be realized with other AFAs. 
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Supplementary material for: Heavy-fermion valence-bond liquids in ultracold atoms: 
cooperation of Kondo effect and geometric frustration 


SU{3) KL MODEL ON A TTL 

In this section we outline a derivation of the KL Hamil¬ 
tonian (3). We start by considering an Anderson impu¬ 
rity model for an isolated triangle in Fig. 1(a), i.e. a 
single-rung case of Eq. (1). In the limit of weak va¬ 
lence fluctuations, we use a Schrieffer-Wolff transforma¬ 
tion to obtain a second-order Kondo-like impurity Hamil¬ 
tonian. Because the periodic Anderson model (1) is addi¬ 
tive w.r.t. triangular rungs, this Kondo impurity model 
can be straightforwardly generalized for a periodic array 
of the “impurities” (triangles), i.e. a TTL. 

The single-triangle Anderson impurity model describes 
three conduction bands coupled to a triangular plaquette 

Haim =Ho + + Ha', (SI) 

Ho = gfcg-Cfc^^Cfcgo-; 
kaa 

+ h-C); 

^ kaa 

Hd =JHSd — CdNd + U ^ n'^nf, 

a^h 

where we used the same notations as in Eq. (1) but omit¬ 
ted the triangle position Xi = i. N is the dimensionless 
system size, e.g. the number of unit cells. 

The last line in Eq. (SI) is the Hamiltonian of an iso¬ 
lated triangle. Due to the strong on-site Coulomb repul¬ 
sion we can treat daa as constrained fermions (no double 
occupancy) with anticommutation relations [1] 

{daaj dbp} = = 0; 

{^aa ; ^5/5} ^ab *^a /3 (f ^a)] ’ 

and consider only three states per vertex a. There are 


hence 4^ = 64 electronic states. The corresponding ener¬ 
gies EdiNd,Sd), 

Edi^, 1/2) = — ed + -^Jh, 

Edi‘2,, 0) = — 2ed + U, 

Ed{2,, 1) = — 2ed + U -\- 2Jdj, 

Ed(3,l/2)=-3ed + 317+^Jff, 

Ed(3,3/2)=-3ed + 317+^Jff, 

are plotted in Fig. SI. For us the most interest has the 
two-electron subspace with a singlet GS. This state is 
three-fold degenerate and corresponds to a global energy 
minimum for U — \ Jh < Cd < 2U -\- ^Jh- The three 
singlet ground states are given in Eq. (2). 

In the strong-coupling regime v Cd, U we can use 
a Schrieffer-Wolff transformation H = e^H aim^~^ [2] 
which is constructed to eliminate iLint in each order in v. 
To the second order we have 

iL2 = i?o + i?d+^[5,Lfi„t]. 

The generator S satisfies the equation [Hd+HQ,S] = 
and has the form [3, 4]: 

^ ^ ■^r^/ka.PddaaPd - h.C. 

Here the summation is over all variable not present in 
the l.h.s. The operators Pd are projectors on the trian¬ 
gle states with energy Ed, and we neglected the electron 
bandwidth in the energy denominators. 

The matrix elements of the commutator in E12 can be 
computed as 


(oi I [5,iLint] 102) = ^ 'y ' _ j^^{0‘l\(dA;aaPddaaPd){d}^^idC-pbll) ~ {d\pCpbp){(^j^^^PddaaPd)\o-2) + h.C. = 


= -y 

N ^ 


I 


1 


Ed{2,Q) - 1 / 2 ) - E(2, 


■ h.C. 


where £1^(2,0) — Ed{i,l/2) = —(2U — ed + aJh) and 
Ed(l,I/2) - Ed{2,Q) = td-U + IJh- We will focus 
on a special case = |C/ when Ed{‘i, 1/2) — Ed{2,Q) = 
Ed(l,I/2)-Ed(2,0) = A. 


Since the states |ai_ 2 ) in Eq. (2) are singlets, the ma¬ 
trix elements (ai| • • • [ 02 ) satisfy the property 

(ni 1 ^ 2 ) — ^a/3 (ni 1 0 . 2 ) — ^a/3(ni|(ia^d^^|ci2) 
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and similarly for {ai\d\pdaa\o, 2 )■ The spin-polarized ele¬ 
ments are given by 


(^1 l^a'I'^5'1' I ^2) — ^ab^a\a2^a‘20. (-> ^ab)') 




SabSaia 2 i^ ^a 2 a ) -I- {l-Sab)SaiaSa 2 b 


Using these expressions and neglecting the constant term 
that arises from Cpb^cl^a = ^abSkpSap - Aaa^pbP^ we ar¬ 
rive at the Kondo impurity model 


^ef — gfeo-Cfc^jg-Cfcacr ^ fl'fa,<^ka' 


kp 


-Yi^ fu 

AT -'a-'' 




kp 


with Vj| = and VI = — from which Eq. (3) is 
recovered by adding a triangle position i, and replacing 
in the last two terms the summation over k and p with 
that over i: 


1 

N 


Yca'a^pacr 

kp 


i 


Note the term gives only a constant energy shift, since 
(ai|ifd|a 2 ) = £1^(2, 0)i5aia2) and therefore was dropped. 
The Hamiltonian Hgf can be reduced to a single-channel 
form using the same arguments as in the main text. 
Hence, in the following we will consider its spin-polarized 
version, omit the spin index a and replace Ska- —t ^k- 

Hel f^kc\^Cka — ^ VabflfbcYcpb (S2) 

k kp 



ea/U 


FIG. SI. Energy levels of an isolated triangle. Nd and Sd 
are the total particle number and spin. In the shaded region 
the GS has Nd = 2 electrons (thick line). Inset: The C^v 
symmetry group of the triangle, which contains the three¬ 
fold axis C 3 and three mirror planes cr«. The numbers 1... 3 
indicate the vertices. 


with Uab = Ul(1 - 6ab) + Vl\Sab- 

Finally we note that an Anderson model similar to 
Eq. (SI) was studied in the context of tunneling through 
triple quantum dots [5, 6], solid-state spin entanglers [7], 
and Kondo physics in frustrated geometries [8, 9], but to 
our knowledge has never been generalized to a lattice. 


PHYSICAL PROPERTIES OF THE SU{3) KL 

In this section we elaborate on the discussion in the 
main text regarding physical properties of the KL model 
(3). First, we will consider the limit of a single SU{3) 
impurity in order to identify parameter regimes relevant 
for Kondo screening. Based on the obtained intuition, 
we then generalize the hybridization mean-field approach 
[10] to handle the SU{3) KL model in Eq. (3). 


Poor man scaling for a single impurity 


We will apply the perturbative (in V|| and Vj_) renor¬ 
malization group (RG), i.e. Anderson’s poor man scaling 
described in Ref. 2 to the Kondo impurity model (S2). 
For simplicity, we assume that the itinerant band has a 
constant density of states po for e G [—D, D]- Notice that 
system is not assumed to be ID. 

The perturbative RG approach is a series of consecu¬ 
tive Schrieffer-Wolff transformations that remove high- 
energy states near the band edges [D — 6D,D] and 
[—D, —D + SD] with SD > 0, SD D. It is assumed 
that there are no electrons (holes) in the first (second) 
interval. We will label states inside (outside) these inter¬ 
vals with momenta q, q' (fc, p), etc. 

The impurity model (S2) can be cast in the form 


Hef — ^kf^kaCka + £ 


qYja^qa 


AT ^abfbfa^kb^pa at ^abfbfa^a'b^qa 

kp q'q 

1 


1 


Y ] Ujf,/j[/a(Cp(,Cgo + CgaCpfc)) 


pq 


where the first two lines contain diagonal terms 
and the last line is an off-diagonal operator Hi = 
- Y Sp 9 yabflfa{cli,Cqa + c\a_Cpb) that transfers electrons 
to and from near the band edges. Note that Hi ~ SD 
because of the q-summation, and we can neglect the term 
with Ylrq'qi’”) ^ The “unperturbed” Hamilto¬ 

nian for the band edges can be taken simply as Hq = 
Trq with C, « ±D. 

The transformation {Hq + Hi)e~^ ~ i£o + 
removes Hi, and the r.h.s. of this expression provides 
a correction to the remaining Hamiltonian which now 
does not contain the band-edge degrees of freedom. All 







3 



FIG. S2. Renormalization group flow deflned by Eq. (S3). 
Solid lines are the trajectories (S4). Labels 1 and 2 correspond 
to the two cases studied in the text, and represent the model 
(S2) with V± = Vj| and V± = — |V|| respectively. 

coupling constants in this new Hamiltonian are defined at 
an energy scale D — SD. The generator S is determined 
by the equation Hi — and can be symbolically 

written as 

E - h.C. 


It is easy to find the RG flow in two special regimes: (1) 
y± = y\\ [full 517(3) symmetry] and (2) V± = — |V||, 
both corresponding to const = 0 in Eq. (S4). Assuming 
that V|| = for some fixed value of 17 = Dg, in the 
first case we have: V||(77) = ~ In ^]. 

In the second case the solution is V||(i7) = V|[°V[l ~ 
In ^]. Upon reducing the bandwidth D < Dq, 

V|| always increases. Hence, for I/||^°^ > 0, the flow 
is always towards strong coupling indicating a Kondo- 
screened GS. The scaling trajectories for a general case 
are shown in Fig. S2. 

The hybridization mean-fleld approach 

The RG flow in Fig. S2 suggests that for > 0, the 
impurity model (S2) always exhibits the Kondo screen¬ 
ing. Similar to the 517(2) case [10], one may expect that 
the same is true for the lattice model (3). Here we demon¬ 
strate this by using a modified slave boson (hybridiza¬ 
tion) mean-field approach that replaces pseudo-fermions 
fia in Eq. (3) with canonical fermions at unit filling: 

i 

where n{ = fj^fia and N is the number of lattice sites. 

The correct slave boson helds can be identihed by di¬ 
agonalizing the interaction part in Eq. (3) on a single 
site assuming that there is one /- and c-fermion. There 
are nine eigenstates 


Here Eq {Eq) are eigenstates of Hq, and Pq (Pq) cor¬ 
responding projectors. The summation extends over all 
relevant variables. 

A straightforward calculation of the commutator 
[5, i7i] yields the RG equations: 


y^{D - SD) =y^{D) + Ui(Ux + 2Ij|) 
y^liD-SD) =Vi^{D) + 3Vl 


or in a differential form: 


dV± 

dlnD 


-Poy±{y± + yi\); 


dyii 

dlnD 


-spoyl 


(S3) 


1"^°^ E //acLivac); (S 6 ) 

IV' 1 . 2 ) =;^(//i4i +‘^'^V/343)|vac); 

IV^3...8) =fL4bW^(^) 

with energies Eq = -(t/|| -|-2Vj_), Ei = 772 = - (V|| - Vj_), 
E 3...8 = 0 , and w = 

The states l'0i,2) transform according to the two- 
dimensional representation of the group C^y [12]. The 
symmetry operations from this group are shown in the 
inset of Fig. SI, and their matrices in the basis {tpip) are 
given by: 


The flow trajectories (integrals of motion) are obtained 
by excluding In 77 from these expressions, i.e. dVi\/dV± = 
3Vx/(2V|| +Ul). This is a Darboux equation [11], whose 
solution is 


(14-ljl)3 



^1 


2 

= const. 


(S4) 


C.= 


uj 0 

0 w* 




0 1 
1 0 




0 w 
w* 0 


ui) = ial) 


The state j'0o) is invariant under all operations from C^y. 

Each of the states (S6) can be associated with a 
Schwinger boson (for a fixed 7), and we can formally 
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write the Kondo term as 

HK=Y.A‘^2'fl,clc^bf^a = 

i 

= '^En{ab'\n){n\ab)fl,cl^,c^bfia = 

in 

= EJnxlnXiri- 

in 

In our case the slave bosons are actually pairing fields: 

XiO ^ ' fiaCia] (S'^) 

^ a 

Xil,2 + W^^/isCis) 

Other Xin with n = 3 ... 8 do not appear due to vanishing 
eigenvalues 

We assume that the system is homogeneous so all av¬ 
erages are i-independent, and perform a Hartree-Fock- 
Bogoliubov decoupling in two “channels”: (i) Hybridiza¬ 
tion channel, by condensing the x-bosons, i.e. consid¬ 
ering the averages Xn = iXm), and (ii) magnetization 
channel, defined by the traceless parts of {fjafib) and 
(claCib)- The latter can be expanded in terms of Gell- 
Mann matrices A* with I = 1.. .8: 

iflaU) + X^ab, 

I 

{4aC^b) = + \5ab- 


with 


Our convention for the Gell-Mann matrices is: 


^0 1 0 
= I 1 0 0 

yO 0 0 

^0 0 1 ' 
A^ = I 0 0 0 
,10 0 



r 

—i 


/ 

II 

i 

0 

o' 

, A3 = 

/ 

^0 

0 

oj 

\ 



0 

-i\ 

/ 

II 

0 

0 

0 

II 

/ 


0 

0 



1 0 ON 
0-10 
0 0 0 ; 

'0 0 0 ^ 

0 0 1 
0 1 Oy 


1 0 0 



yO 0 - 2 y 

The mean-field Hamiltonian has the form: 

Hm¥ = X] [(^fc “ IJ‘c)bab - c\^Ckb + 

k 

k 

~ E [^“'"fcoZ-fca + /-fcaCfea] = 


= E(4a 


'Hk = 


_ f fbc)^ab ^abab 


-V:Sab -^lf6ab+M^^ 


Here summations over color indices are implicit and the 
fc-summation extends over the entire Brillouin zone. The 
chemical potential fXf enforces the constraint (S5), while 
/ic fixes an average number of the itinerant fermions: 

1 




ka 


Finally, the mean helds are = [V||(5ah + Vl(l — 

(5ab)] and V^Va = (V|| -f 2Vj_)xo + (V|| - 

hj_)[xiw + X 2 W*]. The Hamiltonian Hmf can be diag¬ 
onalized by a Bogoliubov transformation [13] 


^ka 

ft 

J —ka 


n+{k) n_(fc) 

= E x%k)xks+ E 


kt 


with n_|_ [n_] is the number of positive [negative] eigenval¬ 
ues Es{k) [—Et{k)] {Es^t > 0) of Hk- These numbers of 
course depend on k and the mean-held parameters. The 
fermionic quasiparticles ■jks and j3kt dehne the particle- 
hole excitations and have only positive energies: 

pn+(fc) n_(fc) 

EE Es{k)xlslks + E 

fc L g j 

We also note a useful relation 
— n-^ = n° — 1 = 


-vE E 


N 


■n+ik) 


7 

4s 


n_(fc) 

- E nL+'n-ik)-8 


Once the quasiparticles are known, we can self- 
consistently determine the order parameters Xn and rn^’k 


Stability of Kondo-screened phases against 
lead-mixing perturbations Hmix 

In the main text we only briehy mentioned effects of 
the lead-mixing term i/mix in Eq. (1) on the heavy- 
fermion phases in Figs. 2 and 3. Here we provide a 
more detailed discussion of this issue. 

The simplest form of i?mix compatible with Csy sym¬ 
metry of the TTL corresponds to hopping of itinerant 
and local fermions around a triangle: 


H ■ — —E ^ 


-A bil^^Ckba Ia kl^^diba + h.C., 

(ab) ,k<T {ab) 


where are hopping amplitudes for itinerant and local 
fermions, and {ab) is a triangle edge (nearest-neighbor 



5 


link). This Hamiltonian is diagonalized by introducing a 
transverse momentum A = 0, 

AiAcr iX(7 

with = -2t^‘^cosA, CfcAa = ^ Ea e“'^“cfca<T and 
similarly for dix^. In all other terms in (1) we 

can simply replace the lead index a with A and re¬ 
peat calculations leading to the KL model (3) using the 
dressed fermions CkXa and diXa- Note, that the replace¬ 
ment a —)■ A also needs to be performed in the ex¬ 
pression (2) for the VB singlets. The d-term in idmix 
is diagonal in the basis of the VB states. For exam¬ 
ple, = (? 2 +? 3 )| 1 )* = “Cill)* because 

Ea^a"* = 0 ’ so in general Y.„d\xad^\<T\a)i = 

Hence, the net effect of idmix amounts to introducing 
Zeeman splittings for c-fermions and 

for the local SU (3) spins: 

5HY^ = Y.{hlc\^c,x + Kflhx). 

iX 

In the context of the SU{2) spin Kondo physics it is 
known [2] that a sufficiently large magnetic field (of the 
order of the Kondo temperature) suppresses the Kondo 
effect. The same will happen in our case: The lead¬ 
mixing hopping amplitudes play the role of an ap¬ 
plied held, and when they become large enough the local 
GS degeneracy on each triangular plaquette will be lifted 


which in turn will destroy the frustrated heavy-fermion 
phases in Figs. 2 and 3. 
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